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Unbroken continuous translational invariance is often taken as a basic assumption in discussions 
of spontaneous symmetry breaking (SSB), which singles out SSB of translational invariance itself 
as an exceptional case. We present a framework which allows us to treat translational invariance on 
the same footing as other symmetries. It is shown that existing theorems on SSB can be straightfor- 
wardly extended to this general case. As a concrete application, we analyze the Nambu-Goldstone 
modes in a (ferromagnetic) supersolid. We prove on the ground of the general theorems that the Bo- 
goliubov mode stemming from a spontaneously broken internal U(l) symmetry and the longitudinal 
phonon due to a crystalline order are distinct physical modes. 
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I. INTRODUCTION 

The low-energy physics of systems with spontaneous 
symmetry breaking (SSB) is dominated by the associated 
soft modes, the Nambu-Goldstone (NG) bosons. The ex- 
istence of at least one such a soft mode is guaranteed by 
the celebrated Goldstone theorem pQ (see also Ref. [5] 
for a comprehensive early review). The question of how 
many NG bosons there are given the symmetry-breaking 
pattern is easy to answer in Lorentz-invariant systems. 
In the general case, it is however difficult; several par- 
tial solutions to the problem have been given over the 
years [SHI] (see also Ref. [6] for a recent review). 

What is common to the general results existing in liter- 
ature is the separate treatment of internal and spacetime 
symmetries. In particular, it is usually assumed that the 
ground state of the system has a continuous translational 
invariance. SSB of translational invariance itself is then 
treated as an exceptional case (an analysis based on ef- 
fective field theory was given in Ref. [?])• Since only a 
discrete translational symmetry is needed to have well- 
defined low-momentum quasiparticles, and thus to dis- 
cuss the NG bosons, such separation seems to be merely 
a matter of convenience. 

In the recent paper [8] , we argued that the distinction 
between internal and spacetime symmetries is somewhat 
artificial, at least as long as SSB is considered. Instead, 
we pointed out that the requirement necessary for stan- 
dard arguments concerning SSB to apply is that the sym- 
metry be uniform. This, roughly speaking, means that 
the infinitesimal symmetry transformation does not de- 
pend explicitly on the spacetime coordinates. Formally, 
it can be encoded in the condition on the translational 
property of the Noether charge density, Eq. ([!]). 
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The objective of the present paper is to complement 
our previous study by a more detailed discussion of spon- 
taneous breaking of translational invariance. Since trans- 
lational invariance is the only uniform spacetime symme- 
try, this completes our analysis and shows that all uni- 
form symmetries can be treated on the same footing. At 
the same time, it provides a generalization of the other 
existing results. What we do assume is an unbroken dis- 
crete translational invariance of the ground state though. 
This is a minimal necessary requirement for the notion of 
NG bosons as soft quasiparticles to be meaningful. One 
can expect that even in the complete absence of transla- 
tional invariance, SSB will give rise to low-lying excited 
states in the spectrum; such a generalization, however, 
goes beyond the scope of this paper. 

We would hasten to emphasize that the issue of spon- 
taneous breaking of translational invariance is, of course, 
not of merely academic interest. Systems where an inho- 
mogeneous state develops spontaneously have been ex- 
tensively discussed in the literature. A classic example 
is the inhomogeneous pairing in superconductors in ex- 
ternal magnetic field, predicted by Larkin, Ovchinnikov, 
Fulde, and Ferrell (LOFF) [S]. Similar type of pairing 
can occur in many other systems, ranging from ultracold 
atomic gases to cold dense quark matter, see Ref. [10] for 
a review and further references. 

In LOFF-type superconductors, the order parameter 
is usually modulated in one dimension. On the contrary, 
complete breaking of continuous translational invariance 
is featured by another prominent example, the super- 
solid [TT]. A supersolid is a state which has both a crys- 
talline order (ordinary long-range order) and a superfluid 
order (off-diagonal long-range order). In the present pa- 
per, we will use it as an example to illustrate the general 
results concerning spontaneous breaking of translational 
invariance. Other examples of systems exhibiting sponta- 
neous breaking of translational invariance include stripe 
ordering in spinor Bose-Einstein condensates [12] and nu- 
merous systems in high-energy physics such as the charge 
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or chiral density wave [13] or the quarkyonic chiral spi- 
ral pij. 

All in all, it is clearly desirable to have a systematic 
framework that can be used to analyze the physics of NG 
modes in systems where translational invariance is bro- 
ken to a (nontrivial) discrete subgroup. Let us remark 
that all our general arguments, presented in Section [TTJ 
are actually valid regardless of whether this breaking is 
spontaneous or explicit. As a byproduct of our inter- 
est in spontaneous breaking of translational invariance, 
we therefore provide a formalism for description of NG 
bosons of spontaneously broken internal uniform symme- 
tries in systems with a discrete periodic structure such 
as crystalline solids or optical lattices [T5] . 

The plan of the paper is as follows. In Section [TTJ we 
discuss in detail the setup and generalize the existing the- 
orems on NG modes, in particular the Goldstone theo- 
rem [T] , the theorem of Schafer et al. [I] , and the Nielsen- 
Chadha theorem [3j. Readers not interested in general 
properties of spontaneous symmetry breaking might pre- 
fer to skip this section and proceed directly to Section |TTTJ 
where we analyze in detail a model for a two-dimensional 
(ferromagnetic) supersolid. As another nontrivial exam- 
ple, we argue how translational properties of the system 
get affected by an external uniform magnetic field and 
discuss phonons in a charged crystal under the magnetic 
field (Section [IV|. In Section [VJ we summarize and con- 
clude. 



II. GENERAL THEOREMS 

As a first step towards the generalization of the ex- 
isting results on NG bosons, we carefully formulate our 
assumptions. This will in particular clarify the scope of 
the validity of our arguments presented in the later parts 
of this section. 



A. Assumptions 

First of all, we should emphasize that when we speak 
of breaking translational invariance, we always have in 
mind spatial translations. Continuous time translation 
invariance, which guarantees energy conservation, will 
be assumed throughout this paper without further men- 
tioning it. As to the spatial translations, we assume dis- 
crete translational invariance in d! directions and con- 
tinuous one in the remaining d — d! directions. This 
requirement is imposed on both the Hamiltonian (or 
Lagrangian) of the theory and its ground state. The 

translation symmetry group is thus x d ' where 



Ci 



(<*') 



{T R = e iP ^ R / h | R = J2ti n * a > ( n i G Z )i and 



{ a i}f=i are the primitive lattice vectors. Furthermore, 
g{d-d') = | e iP 2 . Ax/ft | Ax € R d-d'y where p ig the mQ _ 

mentum operator and Pi , P2 its components in M. d and 



R d d , respectively. Likewise, the spatial coordinate vec- 
tor x can be decomposed as x = (x\, X2) where X\ £ M. d 
and x 2 <G M. d ~ d . Moreover, the vector Xi can be uniquely 
represented as Xi = X\ + R where R is a lattice vector 
and X\ belongs to the Wigner-Seitz unit cell in M. d . We 
will sometimes for the sake of brevity refer to the set of 
all spacetime coordinates, (t, x), using the corresponding 
regular-weight letter, x. 

Apart from translational invariance, the system may 
possess other, internal symmetries at the Lagrangian 
level, which can be spontaneously broken. Following 
the argument presented in Ref. 8j, we always assume 
that such symmetries are uniform, that is, their Noether 
charge density j°(x) satisfies the following condition [16] . 

f(t,x 1 ,x 2 ) 

= e^-^^^T^fiO, XL 0) T R e-^t-iV-aJ/n (1) 

We also assume that the charge associated with a spon- 
taneously broken symmetry (whether internal or trans- 
lational) is time-independent and is given by an integral 
of a local charge density, Q = J d d xj°(t, x). 

Finally, we assume that there is a complete set 
{|£, fe)} of simultaneous eigenstates of the Hamiltonian 
(H\£,k) — ££(k)|£, fc)), discrete translations in M d 
(2k|£, k) — e tkl ' R \£, k)), and the remaining momentum 
operators (P2\£,,k) — Hk2\£,k)). The normalization of 
the states is chosen as (£, fc|£', k') = (2n) d 5 iil S d (k - k f ) 
so that the completeness of the basis is expressed by the 
relation 
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5 |£,fe)<£,fc| = l. 



(2) 



The integration over fci is performed within the first Bril- 
louin zone (FBZ), that is, the Wigner-Seitz cell of the 
reciprocal space^, while the k 2 integration is done over 
the whole M. d ~ d space. The index £ may include dis- 
crete as well as continuous labels such as the band in- 
dex, which together with the momentum uniquely iden- 
tify each eigenstate. 

It is customary to characterize SSB by the number 
of spontaneously broken generators Q a , denoted here as 
7T-BS- This can be usually determined as the difference of 
the dimensions of the symmetry groups of the action and 
the ground state of the theory. However, for the sake 
of general arguments, it is more convenient to use the 
following formal definition: ubs symmetries are said to 
be spontaneously broken if there is as set of (quasi)local 
operators {cf> a (x)}2=\ such that the matrix M given by 
iM ab = (O|[Q O ,0&(O)]|O) is nonsingular. 



B. Goldstone theorem 

The Goldstone theorem 1] claims that whenever a uni- 
form symmetry is spontaneously broken, there should be 
at least one massless bosonic state in the spectrum of 
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the theory. The proof proceeds essentially by finding the 
spectral representation of the commutator of the broken 
charge density j~ and the interpolating field 4> b . Here we 
adapt the standard argument to account for the possi- 
bility of discrete translational invariance and obtain the 
following spectral representation, 

iS ab (uj,k) = J d d +W("*>(0|[j Q V),M0)]|0) 
=2n J2 [% ~ fo(*0) (0\Ta ik ) le, k) (£, fc|0 b (O) |0> 

- s(u + ^ e (-fe))(0|^(0)|e, -*)<£, -fc|ja (fc) |o)] , 

(3) 

which holds provided fci lies in FBZ. This expression is 
formally identical to one assuming full continuous trans- 
lational invariance; all effects of only discrete transla- 
tional symmetry in M. d are hidden in the definition of 
the charge density averaged over the unit cell, 

Ta (k) ^^-f j o (0,xx,0)e- ife ^, (4) 

"uc Juc 

where f2 uc is the volume of the Wigner-Seitz unit cell, 
over which the X\ integral is performed (this is indicated 
by the subscript "uc"). Note that in the limit of f2 uc — > 0, 
ia(fe) Ja(0) which leads to the usual spectral represen- 
tation with the full continuous translational invariance. 

The proof of Eq. ^ is straightforward but somewhat 
tedious, and we thus refer the reader to Appendix [A] for 
details. Consider now the limit |fc| — > of the commuta- 
tor Fourier transformed only in the spatial coordinates, 

Urn iS ab (t,k) = lim [ d d xe- lk - x (0\[j°(x),M0W) 
|fc|->o |fe|->o7 

= <0|[Q«,&(0)]|0)=iM o6 . (5) 

Using the assumption that the charges Q a are time- 
independent (which is most naturally ensured when 
j'2 satisfies a continuity equation together with a spa- 
tial current j a , and the corresponding surface term in 
the integration over space vanishes), we deduce that 

lim iS ab (cj,k) cx 5(ui). Combining this with the spec- 
|fcH-o 

tral representation ([3]), we find 

M ab = 2 Im[(O|^,O)(e,O|0 b (O)|O)] , (6) 

£|e 4 (O)=0 

where j a = j a{0) . 

The condition of SSB (detM ^ 0) ensures that the 
right-hand side of this equation is nonzero for some a, b. 
One might at first think that this is due to the contribu- 
tion of degenerate ground states obtained by the action 
of broken symmetry transformations on |0). However, 
this contribution drops at large volume £1 asymptoti- 
cally as 1/Q since such degenerate ground states only 
exist for isolated values of momentum (for which k2 = 



and ki is equal to some of the vectors of the reciprocal 
lattice) [21 [17]. It then follows that there must be at 
least one intermediate excited state \n,k) such that: (i) 
(0|j' a ( fe )|n, k) and (n, k\4> b (0)\0) ^ for some a, b in 
the vicinity of k = 0, and (ii) e n (0) = 0. Such inter- 
mediate states are exactly the Nambu-Goldstone modes 
associated with spontaneous breaking of the charges Q a . 
Denoting \n) = \n, 0) and the number of these states as 
rtNG, Eq. ^ becomes 

"NG 

M ab = 2j2 Im[(O|j^n)(n|06(O)|O)]. (7) 

n=l 

Before concluding the discussion of the Goldstone the- 
orem, let us remark that in case there are long-range 
interactions in the system, the surface term appearing 
in the space integral of the continuity equation for the 
Noether charge density and its current need not vanish. 
The integral charge then depends on time and massive 
modes can contribute to Eq. ([6]). This is why the longi- 
tudinal mode in a Coulomb-interacting Wigner solid (in 
three spatial dimensions) acquires an energy gap, equal 
to the plasma frequency [18]. 

Finally, note that the above presented proof of the 
Goldstone theorem relies essentially on the operator for- 
malism and the translational property ([T]) of the charge 
density as well as the time independence of the integral 
charge. In concrete applications, it is often advantageous 
to use an alternative proof based on the effective action 
formalism ^ I19j , which provides a direct connection to 
the symmetries of the corresponding classical Lagrangian 
system. For the sake of completeness, we sketch a modifi- 
cation of this proof suitable for systems with only discrete 
translational invariance in Appendix [B] 

C. Theorem of Schafer et al. 

In the preceding subsection, we carefully formulated 
the Goldstone theorem as showing the existence of at 
least one NG boson. The natural question to ask is, of 
course, whether anything can be said in general about 
their actual number, tt-ng- A common lore in textbooks 
on relativistic field theory is that assuming Lorentz in- 
variance, rtNG equals the number of spontaneously bro- 
ken generators, jibs- Nevertheless, there are numer- 
ous examples of Lorcntz-noninvariant (both relativistic 
many-body as well as intrinsically nonrelativistic) sys- 
tems where this simple relation does not hold (see Ref. [6] 
and references therein). 

We start our discussion with the result of Schafer et 
al. [4] who formulated the following theorem [20 : Pro- 
vided that (0|[<2 a , <2i,]|0) = for all a, b, the number of 
NG modes ting is equal to the number of broken symme- 
tries riBS- Here we generalize this theorem by taking into 
account the possibility of a discrete translational symme- 
try in M. d (whether as a result of spontaneous breaking 
of continuous translational invariance or of the intrinsic 
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setup of the system). We follow essentially the argument 
of the review [6] . 

Let us first show that ti^g > "bs- Based on the re- 
maining translational symmetry of the ground state, we 
infer that (0|[Q a , Q&]|0) oc (0| [Q a , jb] |0) and, additionally, 
using Eq. Q, 
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Q a ,jb]\0) =2i Vim [<0|j»<n|j 6 |0>] 



(8) 



Defining an nes x "ng matrix A as A an — (0\j a \n), the 
assumption that (0|[Q a , Qh]|0) = for all a, b together 
with Eq. Q implies that AA< is Hermitian as well as real. 
It can thus be diagonalized by a real orthogonal matrix 
O such that A 1 A A = OAA^O T = diag(Ai, . . . , A„ BS ), as- 
suming without lack of generality < Aj. < ■ • • < A„ BS . 
Now if the number of NG bosons were smaller than nes > 
the rank of A, and hence also the rank of AA\ be- 
ing at most riNGj would have to be smaller than nes 
as well. The smallest eigenvalue Ai would then neces- 
sarily vanish. However, since the orthogonal transfor- 
mation by O can be interpreted as a change of basis of 
the broken generators, Q' a — O a bQb, this would lead to 
(0|ji |n) — (n\ji |0) = for all n, in contradiction with 
the broken symmetry assumption, det M ^ 0. (Note that 
A f = OM so that determinants of M and M' are equal 
up to a sign.) Thus, the inequality ting > ^BS follows. 

In order to prove the opposite inequality, tt-ng < ^bs, 
observe that A^A is positive semidefinite and Hermitian, 
and can therefore be diagonalized by a unitary matrix 
U as A'^A" = WA^AU = diag(pi, . . . , p„ NG ), again as- 
suming without lack of generality that < p± < ■ ■ ■ < 
p„ NG . The transformation A —> A" = AU corresponds to 
a redefinition of the basis of states, \m") = \n)U nm , since 
A'am = A an U nrn = {Q\j a \m") (summations over repeated 
indices are implied). Now if the number of NG bosons 
were larger than jibs, we would analogously deduce that 
(0\ja\ m ") = for all a and for m = 1,2, ... , ting — "bs- 
However, these zero modes do not couple to any broken 
currents and thus do not satisfy the very definition of a 
NG boson. 

Putting both pieces together, we have proven that the 
number of NG modes is equal to the number of broken 
symmetries if (0|[Q a , Q&]|0) = for all a,b. 



D. Nielsen— Chadha theorem 

Nielsen and Chadha [3] discovered that the number 
of NG bosons is tightly related to their dispersion rela- 
tions. Concretely, they showed that under certain tech- 
nical assumptions, the dispersion relation of a NG boson 
is proportional to a positive integer power of momen- 
tum in the long- wavelength limit. Classifying NG bosons 
with an odd (even) power of momentum in the dispersion 
relation as type-I (type-II) and denoting their numbers 
as 7ii and ran, respectively, they proved the inequality 
n\ + 2n\i > ubs- In order to generalize their result to 



systems with only discrete translational symmetry, we 
basically just need to replace the spectral representation 
in their original paper by Eq. pi). We shall nevertheless 
present some details of the argument. 

Let us define the vectors (A n )b — (0\<j>b(0)\n) and v a = 
A n (n\j a \0}- Let in addition p be the number of linearly 
independent vectors among the v a 's (obviously p < ting)> 
and set An = nes ~ P- By construction, there is a set of 
coefficients C" (a = 1, . . . , An) such that G^v a = and 
simultaneously rankC = An. Using Eq. ([3]), we obtain 

iC^S ab (LJ,k)^ 

27rV/(o;- i £n (fe))<0|C^ (fe) |n,fc)(n,fe|^(0)|0) (9) 



in the vicinity of k = 0. Here we dropped the contribu- 
tion from massive modes and the second term of Eq. (|3| , 
since they both vanish at k = 0. Assuming along with 
Nielsen and Chadha the exponential decay of (expecta- 
tion values of) commutators of local operators separated 
by a large spacelike interval, we infer that S a b(t,k) is 
an analytic function of k. In other words, the support 
of S a b(oj, k) should be differentiable with respect to k. 
As a consequence, for any \n,k) that appears in Eq. ^ 
for some a and b, s n (k) must be an analytic function of 
k. We will call such NG modes type-C and denote their 
number as nc- 

As the next step, we observe that the An vectors 
C"*v a = A n (n\C2*j a \0) must be linearly independent. 
Indeed, if this were not true, there would be a set of coeffi- 
cients /3 a such that /3 a C"*i> a =0 and thus /3 a C™*M ab = 0, 
implying that M has a zero eigenvalue, in contradiction 
with the assumption det M ^ 0. From the fact that the 
vectors C^*v a are linearly independent, we can imme- 
diately conclude that the number of NG states |n) for 
which (Q\C"j a \n) ^ for some a, must be at least An, 
hence nc > An. 

The remaining NG modes, for which (0|C"j a |n) = 
for all a, will be called type-N and their number denoted 
accordingly as %. The above argument does not con- 
strain in any way the dispersion relation of type-N NG 
bosons; in particular, it is allowed to be linear, e(fe) oc |fe|. 
Using the fact that, by construction, + nc = jingj we 
obtain the inequality 

n N + 2n c > n NG + An = n B s + ("ng ~ p) > ™bs- (10) 

Note that in Ref. [S], we actually proved the equality, 
tin + 2nc = nes, using a slightly different definition of 
the two types of the NG modes. 

In fact, the derivation of the inequality ( 10 1 does not 



make any use of rotational invariance which was implic- 
itly, and somewhat unnecessarily, assumed in Ref. [3]. 
Once this constraint is released, the classification of NG 
bosons can no longer be based on a mere power of the 
momentum |fc| in the dispersion relation. A convenient 
generalization is to define as type-II those NG modes 
whose dispersion is analytic around k = 0, and as type- 
I the remaining ones; this classification reduces to that 
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of Nielsen and Chadha once rotational invariance is re- 
stored. We have thus shown in other words that every 
type-C NG boson is simultaneously type-II. Using the in- 
equality ( 10 ), we then find n\ + 2n\\ > tin + 2uq > "BSj 



which reproduces the Nielsen-Chadha inequality. 



E. Implications for spontaneous breaking of 
translational invariance 

In the preceding subsections, we saw how the general 
theorems known in literature can be extended to account 
for the possibility that the ground state of the system 
has only a discrete translational symmetry. We are par- 
ticularly interested in the situation when this arises as a 
result of a spontaneous breakdown of a continuous trans- 
lational invariance. We will therefore devote this subsec- 
tion to some discussion of systems with the symmetry 
breaking pattern 



Sint X Q%> -> G'u 



x G (d ) x G {d - d,) 

int x »D x »C ' 



(11) 



where Q m t describes the internal symmetry group, which 
can also be (partially) spontaneously broken. 

First, note that the Noether charge of translational 
invariance, that is, the momentum operator, commutes 
with all other broken charges Q a . If we for a mo- 
ment assume, as in the theorem of Schafer et ai, that 
(0| [Q a , Qb] |0) = for all pairs of broken internal symme- 
try generators, we immediately conclude that the number 
of NG bosons equals the number of broken generators. 
Therefore, there is one NG boson for each spontaneously 
broken translation generator as well as for each sponta- 
neously broken generator of an internal symmetry. 

In the general case, it is convenient to define the an- 
tisymmetric matrix p by ip a b = lhn i7 _1 (0|[(5 a , Qfc]|0) 

where a, b = 1, . . . , nes- In the recent paper [H], we pro- 
posed that the number of NG bosons is generally related 
to the rank of the matrix p by 



?ibs - «ng = \ rank p. 



(12) 



The fact that translation generators commute with ev- 
erything else now implies that rank p — rank p lnt , where 
Pint is the analogous matrix defined using the generators 
of broken internal symmetries only. We can thus infer 
that the counting of NG bosons of spontaneously broken 
translational invariance is not affected by the presence of 
other spontaneously broken, internal symmetries: each 
spontaneously broken translation generator gives rise to 
exactly one NG boson. In Section |III[ we will discuss 
examples of both types of systems where rank p is zero 
and nonzero. 

So far, we have tacitly assumed that the generators 
of space translations commute with one another as well 
as with other symmetry generators. Nevertheless, in 
presence of external fields, this property can be violated 
while still maintaining well-defined translational symme- 
try. Let us consider, for instance, a system of (inter- 
acting) charged particles in a uniform external magnetic 



field B. Due to the explicit dependence of the vector 
potential on the coordinates, the action is no longer in- 
variant under the usual translations. Instead, it is in- 
variant under magnetic translations, that is, translations 
combined with a gauge transformation of the electromag- 
netic field. The momentum operator P as a symmetry 
generator is then replaced by the operator P B , which 
satisfies 

[if, if] = -ihqe ljk B k Q , [Pf, Q ] = 0, (13) 

where q is the charge of the particles and Qo the U(l) 
Noether charge, that is, their total number. Note 
that the rank of the matrix p B defined by ip B = 

lim FT 1 (0\[Pf , P?]\0) equals two if the magnetic field 

as well as the expectation value of Qo is nonzero (other- 
wise the rank is zero). Eq. ( 12 ) then suggests that in case 
magnetic translational invariance is spontaneously bro- 
ken, there is one less NG boson than naively expected 
based on the number of broken translation generators. 
We will see an example of this in Section IV [2"T] . 



III. NG MODES IN A SUPERSOLID 

In this section, we will illustrate the general results 
obtained above on the example of NG modes in a (fer- 
romagnetic) supersolid in d spatial dimensions. We will 
employ the mean-field approximation at zero tempera- 
ture. 



A. Formalism 

1. Definition of the model and its symmetry 

Let us consider the class of theories, defined by the 
following Lagrangian, 



= J d d x^\{x) 
-\ f d d xd d y^\{x)^i{x)V{x - y^MVjiv), 



h 2 v 2 

2m 



(14) 



where 'J = (^i, . . . , ^n) is an iV-component complex 
scalar field, and the potential V(x) — V(—x) repre- 
sents a finite-range interaction. The symmetry of the 
Lagrangian is obviously U(iV) x Gq • The simplest case 
N = 1 corresponds to the usual Bose-Einstein conden- 
sate (BEC) without internal degrees of freedom, while 
the widely discussed N = 2 case can be understood as a 
spin-1/2 BEC (see, for instance, Ref. 22 ). 

It is convenient to analyze the model using dimension- 
less variables; this makes the results universally applica- 
ble to different physical systems, as long as they share 
the same set of symmetries and degrees of freedom. We 
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first introduce a characteristic length scale of the poten- 
tial, a, in terms of which the strength of the interac- 
tion can be measured by the parameter Vq, defined by 
a d Vo = J d xV(x). Next, we rescale the spacetime co- 
ordinates as t — > thjii and x — > xa. (For conventional 
reasons we use the same symbols for the new, dimension- 
less coordinates; as we do so systematically from now 
on, no confusion can arise.) Finally, we introduce the 
dimensionless potential, v(x) — V(x)/Vq, and the field 
iJj(x) = 'i!(x)y // Voa d / [i. The dimensionless Lagrangian, 
L[ijj\ = L[^]Vo//i 2 , thus becomes 



L[i/>] = J d d x4{x) \id t + + <Pi( x ) 

■i J d A xd^yi}\{x)i)i{x)v{x - y)ip](y)ipj(y), 



(15) 



where g is a dimensionless coupling parameter defined as 
g = p,(h 2 / ma 2 )^ 1 . 



2. Classical field 

The first step in the mean-field analysis is to find the 
ground state, that is, the static classical field configura- 
tion ipo( x ) which minimizes the energy, 



E[tp] = J d d x 



1 

2f) 



Vi()*(x) ■ Vtpi(x) - tf(x)tl)i(x) 



+ - J d d xd d yil)*{x)il) l {x)v{x-y)ip*{y)ip :j {y). 

(16) 

To that end, we exploit the symmetry of the problem and 
cast the spinor tpo(x) in a "canonical" form. One should 
nevertheless be careful since the field in general depends 
on the coordinates. Note that it can always be written 
as ip(x) — U(x)ip'(x), where U(x) is a unitary matrix 
and ip'i(x) — Sn4>(x) with real 4>{x). As long as the field 
ip(x) does not contain any vortices, both U(x) and V>(a:) 
are differentiable. It then follows that 



E[ip] - E[ijj'] = 



N 

E 



, rP(x) 2 \VU a (x)\ 2 



d d x 



2g 



> 0. (17) 



The equality occurs only if Un is coordinate-independent. 
Thus, by a global U(JV) rotation, ipo(x) can always be 
chosen as ipoi(x) = Sn^^x) with a real field $o{x). Field 
configurations with vortices tend to have a higher energy 
and we therefore assume we can discard this possibility. 

One (nontrivial) stationary point of the energy func- 
tional ( 16 ) is always given by a uniform field. Due to 
our normalization, J d d xv(x) = 1, this corresponds to 
ijo(x) = 1 with the energy density —1/2. The uniform 
solution indeed minimizes energy for small g. However, 
when g becomes large, the energy cost of creating spatial 
modulation of il>{x) is suppressed. Depending on the de- 
tailed structure of the potential v(x), this cost can then 



be overwhelmed by the energy gain from the nonlocal 
interaction. As a result, inhomogeneous field configura- 
tions may become energetically favorable. 



3. Spontaneous symmetry breaking 

As soon as the classical field tpo( x ) assumes a nonzero 
value, the symmetry of the Lagrangian is spontaneously 
broken. Owing to the fact that the spinor ipo( x ) points 
in the same direction everywhere in space, tpoi(x) — 
Snipo(x), the internal U(iV) symmetry is broken down to 
its U(7V — 1) subgroup; 2N — 1 of its generators are thus 
broken. Therefore, when ^)q{x) is homogeneous the resid- 
ual symmetry is U(7V — 1) x and n B s = 27V — 1. We 
will call this solution the superfluid phase. On the other 
hand, if ipo( x ) is inhomogeneous and periodically modu- 
lated in all spatial directions (the supersolid phase), the 
symmetry of the ground state is merely U(7V — 1) x 
and n B s = 2N - 1 + d. 

To predict the number of NG modes, we have to 
evaluate the rank of the matrix p in Eq. (12). To 
that end, note that the 2A^ — 1 broken internal sym- 
metry generators fall into three classes: the generator 
of phase transformations of the first component of the 
spinor, (Ti)jk — <5ji<Jfci/2; the N — 1 real generators 
(Tf-) jk = (SkiSje + 5jiS M )/2; the N - 1 pure imag- 
inary generators {Tj)jk — i(Skifij£ — o~ji8ki)/2, where 
£ = 2, 3, . . . , N. Using these matrices, the broken gen- 
erators can be represented on the Hilbert space of the 
system by the operators Q a — J d d xipl(x)(T a )iji/jj(x). 
Within the mean-field approximation, the fluctuations 
are neglected and one easily derives 



n 



(o\[Qf,Ql]\o) 



iSft 



2il, 



(0|Qi|0) 

d d xij {x) 2 



(18) 



for I' = 2, . . . , N and (0| [Q a , Q b ] |0) = for other com- 
binations. We thus conclude that (l/2)rankp = N — 1. 

In the case N = 1 (usual BEC), we can now use the 
theorem of Schafer et al. to assert that riNG — "bs 
1 in the superfluid phase and %g = n^s = 1 + d in 
the supersolid phase. For N > 2, Eq. (12) implies that 
"ng = N in the superfluid phase and 7Jng — N + d in 
the supersolid phase. To be consistent with the Nielsen- 
Chadha theorem, at least AT — 1 of the NG modes must be 
type-II. This is in agreement with the result of Ref. [53] 
which analyzed the relativistic version of the model with 
a local interaction, V(x) oc S d (x). The detailed nature 
of the NG modes in our model will be investigated in the 
following subsection. 
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4- Elementary excitations 

In the mean-field approximation, the excitation spec- 
trum is found simply by expanding the Lagrangian to sec- 
ond order in fluctuations around the classical field. We 
parameterize the field as ip(x) = ipo( x ) + <P(x) + itp{x) 
and plug this back into Eq. (15). Up to second or- 
der in the fluctuation fields fa ip, the Lagrangian reads 
L[ip] = L[ip ] + L 2 [fa <p;ipo]-\ , where 



i[(f>,ip',ipo] = j d d x ^fa{x)JC x fa(x) + ipi{x)K. x ipi{x) 
+ [ipi(x)d t (j>i(x) - 4>i{x)d t (pi(x)] | 
d d x d d y fa(x)w(x,y)fa(y) 

(19) 

up to a term which is a total time derivative. In the above 
equation, JC X — (V 2 /2g) + 1 — u(x) is the kinetic term, 
u(x) = J d d yv(x — y)\i>o{y)\ 2 1S the effective periodic 
mean-field potential due to the background ipo(x), and 
w(x,y) = ipo(x)v(x — y)'tpo(y) captures the effects of the 
non-local interaction. 

Physically, the fluctuations {</>i}fl 2 an< ^ {fi}iLi corre- 
spond to NG fields of the spontaneously broken internal 
symmetry. Namely, the {fa}^L 2 are obtained from the 
classical field tjjo( x ) by real rotations, generated by the 
matrices T?, while the {<Pi}fLi are induced by symmetric 
unitary transformations, generated by the matrices Xf 1 
(or T\ in case of i = 1). On the other hand, the fa(x) 
mode corresponds to displacements of the classical field 
configuration, ijjo{x + — V-'o( £C ) — €( x ) ' V^oC 33 )) as 
well as to density fluctuations. This can be seen by not- 
ing the role of <f>i{x) in the density-density correlation 
function, 



(Q\T{6n(x)6n(y)}\0) 

~ 4Mx)MymT{Mx)My)M 



(20) 



where Sn(x) = ipj(x)ipi(x) — (0\ipj (x)ipi(x)\0i) , and the 
ellipsis stands for terms of higher order in the fluctuation 
fields. 

The dispersion relations of the various excitation 
branches are determined formally by a diagonalization of 
the Lagrangian ( 19 \. In practice, we accomplish this task 



by solving the Euler-Lagrange equations for the fluctua- 
tion fields, 



(21) 



d t ifi(x) = K x fa(x) ~ 2Sn / d d yw(x,y)4> 1 {y) 
d t (j>i{x) = -JC x ipi{x). 



Note that the equations of motion for fa and tpi are cou- 
pled; we will see the importance of this coupling later. 
Based on the residual discrete symmetry, we can 

assume the normal modes of fa,tpi to have the Bloch 
form, that is, 



i , n (k,x)=J2fa, n}G (k)e i ( k + G >- 



,(fc)t 



(22) 



(analogously for ipi) with a band index n and a crystal 
momentum k S FBZ. The solution Lo n (k) represents the 
allowed excitation spectrum. In particular, those modes 
which satisfy w„(0) = are NG modes, assuming they 
couple to the spontaneously broken currents. 



B. Numerical results 

As a concrete illustration of the general setup, we 
choose the simplest model for a supersolid, defined by 
the the step-function potential V(x) = (Vo/7r)0(a — |jc|) 
in two spatial dimensions |24j . Below, we summarize the 
results of our numerical computations. In order not to 
obscure the physics with technicalities, we describe some 
details of the methods we used in Appendix [C] In partic- 
ular, in Appendix |C 1| we explain ho w th e classical field 
^oix) was found, and in Appendix |C 2| how the band 
structures presented in Figs. 2(a) and 2(b) were obtained. 



1. Superfluid phase 

For g < go rs 38.4, the energy ( |16[ ) is minimized by the 
homogeneous solution ^q{x) — 1 25 . In this superfluid 
phase, the Schrodinger equations ( plj ) are solved straight- 
forwardly by going to the energy-momentum space, tak- 
ing the form Ai(k) (fa(k),(fi(k)) = (separately for 
i = 1, . . . , N), where 



Ai(k) 



{k 2 /2g) + 2S llV {k) 
iuj(k) 



—iui(k)\ 
k 2 /2g ) 



(23) 



Here fa(k) = J d 2 x fa(0,x) e - lk - x+lul ^ t [analogously 
for <pi(k)] and v(k) = / d 2 x v{x) e - ik x = 2J 1 (|fe|)/|fe| 
[J n (x) is the Bessel function]. From the condition 
det Ai(k) = we find the quasiparticle dispersion re- 
lation, uji(k) = y/(k 2 /2g) [(k 2 /2g) + 2S a v(k)\. Around 
k = 0, the dispersions and eigenvectors corresponding to 
the individual modes are expanded as 



wi(fc) ~ \k\/y/g, (^,^)~(0,1) 



(24) 



io 3 (k)~k 2 /2g, 



G 



to leading order in momentum. We illustrate these dis- 
persion relations in Fig. 1(a) for the coupling set to the 
transition point, g = go — 0. In summary, there is one 
type-I Bogoliubov mode and N — 1 typc-II modes with 
a quadratic dispersion and a circular polarization, which 
can be interpreted as ferromagnetic magnons (precession 
modes) (26) . This conclusion is in exact agreement with 
the prediction we made in Section |III A3| 

It is worth emphasizing that while at k = the Bogoli- 
ubov mode corresponds to a pure phase fluctuation, tp\, 
at nonzero momentum it is given by a mixture of tp± and 
fa. In other words, the Bogoliubov mode generates not 
only a phase modulation but also a density modulation of 
the condensate "00- This is why we can associate it with a 




FIG. 1. (color online), (a) The dispersion relations of NG 
bosons in the superfluid phase at g = go — 0. The linear 
branch is the usual Bogoliubov mode, while the quadratic 
branch (with N — 1 fold degeneracy) can be interpreted as 
a ferromagnetic magnon. The dispersion of the Bogoliubov 
mode has a roton-like minimum around |fc| ~ 4.80, which 
corresponds to the magnitude of the reciprocal lattice vector 
at g — go + 0. (b) Three-dimensional plot of the classical 
field ipo(x) in the supersolid phase. The energy functional is 
minimized by a hexagonal (triangular) lattice with the lattice 
constant approximately 1.50 at g = go + 0. Note that tpo(x) 
does not have any nodes. 



density wave and detect it in Bragg spectroscopy experi- 
ments |57j . However, this very fact has also caused some 
confusion regarding the independence of the Bogoliubov 
mode and the longitudinal phonon originating from the 
crystalline order in the supersolid phase. We will shortly 
see that these two are, indeed, distinct modes. 

Finally, let us remark that the above discussed mix- 
ing can be traced back to the term ^d t ip in the La- 
grangian ( 15 ). A similar term also arises in the analogous 



relativistic model where it is induced by the chemical po- 
tential |23) . On the other hand, in a Lorentz- invariant 
model where SSB is triggered by a "wrong sign" of the 
mass term of the field, no such mixing would arise, and 
the 4>i and ifi modes would be completely decoupled. 



2. Supersolid phase 

For g > go, periodic spatial modulation of the order 
parameter 4>q{x) is energetically favored. The transi- 
tion at g = g is of first order within the mean-field 
approximation [24] , and can be roughly understood as 
Bose-Einstein condensation of field modes with momen- 
tum determined by the reciprocal lattice vectors. This 
interpretation is supported by the roton-like minimum in 
the dispersion relation of the Bogoliubov mode shown in 
Numerical computation, following the strategy 



Fig. 1(a) 



outlined in Appendix |C 1[ shows that the solution with 
the lowest energy has the form of a hexagonal lattice, 
plotted in Fig. |l(b)| 



By solving the coupled Schrodinger equations (21 1 us- 
ing the method explained in Appendix |C 2| we found the 
band structures in the supersolid phase at g = go + 0, 
shown in Figs. 2(a) and 2(b) This result was obtained 
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FIG. 2. The numerical result for dispersion relations in the su- 
persolid phase at g = go + 0- Panel (a) displays the excitation 
branches in the <f>i, (pi sector. We can see three linear disper- 
sions around the F point, which correspond to the Bogoliubov 
mode and two phonons. The inset shows the irreducible part 
of FBZ. Panel (b) displays the excitation branches in the <f>i, 
ifii (i > 2) sector. The dispersion relation of the gapless mode 
is quadratic at low momentum, as can be seen in the inset 
which zooms in the vertical axis. 



using the first 73 G vectors around the origin of the recip- 
rocal lattice. The convergence of the result was checked 
by a comparison with one using only 61 G vectors. 

In the i = l sector, three of the excitation branches 
are linear around the T point (the origin of the fe space) , 
corresponding to a longitudinal phonon, a transverse 
phonon, and a Bogoliubov mode. On the other hand, 
each i > 2 sector contains just one massless mode with a 
quadratic dispersion relation — the magnon. In summary, 
there are in total N + d NG modes, N — 1 of which are 
type-II, exactly as predicted in Section [ill A 3| 

In addition, we have checked numerically that the dis- 
persion relations of all the gapless modes are isotropic 
at the r point; the NG modes of the model will thus be 
described by a rotationally invariant low-energy effective 
field theory. Our findings are consistent with previous 
works for the case N = 1, based on the effective La- 
grangian method |28) . A low energy effective field theory 
describing simultaneously the Bogoliubov mode and the 
phonons has recently also been worked out, in a different 
context, in Ref. [2"5] , 

We would like to emphasize that while the specific re- 
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suits presented in Figs. 2(a) and |2(b)j rely on the mean- 
field approximation, the qualitative conclusion that the 
Bogoliubov mode and the longitudinal phonon are dis- 
tinct physical excitations is exact. Following the gen- 



eral argument of Section II E in the case of N = 1 it is 
based on the theorem of Schafer et al. and thus can be 
considered rigorously proven. For N > 2, it is strongly 



supported by the conjectured equality ( 12 ) 



IV. 



PHONONS IN A CHARGED CRYSTAL IN 
EXTERNAL MAGNETIC FIELD 



As anticipated in the discussion in Section HE the 



spectrum of NG modes can become more complicated 
once the generators of space translations do not commute 
with each other as well as with generators of other, in- 
ternal symmetries. To illustrate this subtlety, we shall in 
this section discuss briefly a simple quantum-mechanical 
model for a charged crystal. 



A. Quantum-mechanical model 

The Hamiltonian for N particles with the common 
charge q in non-relativistic quantum mechanics is given 
by 



H = 



N 

E 



\Pi 



qA(xj 
2m 



~J2 v ( x i~ x j)> 



(25) 



where A(xi) is a vector potential which determines the 
external magnetic field B(x) via V x A(x) = B(x), and 
v(x) is a repulsive interaction potential which is assumed 
to be invariant under space inversion, v(x) = v(—x). 
This model, with the Coulomb potential v(x) = e 2 /\x\, 
has been used in the context of the Wigner solid [50] , 
Nevertheless, since our goal here is to discuss the conse- 
quences of SSB in the presence of a magnetic field, we 
will assume that the potential v(x) has a finite range. 

The canonical coordinate and momentum Xi a ,Pjp 
are required to satisfy the usual commutation relation. 
(Within this section, we use Greek letters to denote 
spatial indices.) Because of the vector potential, the 
operator of total momentum J^ili Pi n0 longer com- 
mutes with the Hamiltonian. When the magnetic field 
is uniform, we can instead define the operator pi = 
Pi - q [A(xi) - B X Xi\. The sum, = £\ =1 p ia , is 
conserved and plays the role of the generator of spatial 
translations ([iP^,Xi^\ = fr5 a p)- At the same time, it 
satisfies the commutation relation (13) with Qq = N. 



lattice node, x, = Ri. The spectrum of oscillations 
of the crystal lattice can be determined by resorting 
to the harmonic approximation, in which the Hamilto- 
nian takes the form H = 5Z i=1 [pi — qA(xi)} 2 / (2m) + 

Xi — Ri and 



(m/2) W°fu ia Ujp, where u { 



mW?f 



-d a dfjV 



X — Ri 



N 



d a dpv 



x — H , , 



(Mi), 
(* = .?)• 



(26) 



(We have also introduced the notation = Ri — Rj.) 
Let us now assume that the total particle number N ex- 
tends to infinity or that periodic boundary conditions are 
imposed. Then the Heisenberg equations of motion — 
unlike the Hamiltonian itself — have the lattice transla- 
tional symmetry, since only the magnetic field B ap- 
pears in them. Thus, by using the Fourier compo- 
nents, = J2iW"/e~ lk ' Rij for any fixed j, and 
Uk a (t) — J2i u la e l " ( - k ' )t ~' lk Ri , we simplify the equations 
to the form M a p(k)ukp — with M a p(k) — 8 a puj(k) 2 — 
ie a /3 1 uj(k)b' y — and 6 7 = qB 1 /m. As before, the con- 
dition det M (k) — determines the dispersion relations 
of the normal modes of the crystal lattice. 

For the sake of brevity, we shall consider the simplest, 
albeit unrealistic, model — the isotropic crystal in three 
spatial dimensions, where 



(^-^)^(|fe|) 2 + ^HW) 2 - (27) 

In the absence of the magnetic field, the functions aj f (|fc|) 
and cj£(|fc|) give the dispersion relations of transverse 
and longitudinal phonons, respectively, and it is there- 
fore natural to assume both of them to be linear at low 
momentum. The determinant of M(k) is easily found, 



detM ={uj z -ujf) z (uj z -uf) -lj z (lj z - uj z t )b z sin' 
-uj 2 (uj 2 -ujhb 2 cos 2 8, 



(28) 



where 9 is the angle between the vectors k and b, and 
we have for simplicity omitted the argument k where 
necessary. At low momentum, the dispersion relations of 
the three normal modes can be evaluated analytically, 



gapped mode, 



sin 2 9 + oj 2 cos 2 6, type-I NG mode, (29) 



LU 2 LOe 



\b\\/uj 2 sin 2 9 + cuj cos 2 i 



z, type-II NG mode. 



B. Phonon spectrum 

Suppose that we have found a stable crystal config- 
uration where the i-th particle is localized around the 



In order to determine the polarization of the modes, that 
is, the orientation of the vector Uk, it is convenient to 
rewrite the equation of motion in a vector form 

oj 2 u - iuju x b - (bJ 2 P t u + uj 2 P e u) = 0, (30) 
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where Pt : i are projectors to transverse/longitudinal di- 
rections with respect to momentum k, defined by the 
respective terms in Eq. (l27|). For the gapped mode, the 



terms in Eq. (30) proportional to uif e can be neglected, 



which leads to the condition u = iuxb/\b\. In words, this 
mode is circularly polarized in the plane perpendicular 
to the magnetic field, in the opposite direction than the 
usual cyclotron motion. For the type-I NG mode, we like- 
wise neglect terms in Eq. ( 30 1 proportional to us 2 as well 



as uif t to arrive at the constraint u x b — 0, which simply 
means that this mode is always polarized along the di- 
rection of the magnetic field. Finally, for the type-II NG 
mode we neglect the term in Eq. ( 30 ) proportional to ui 2 
and obtain the condition uxb = (J/uj)(uj 2 PtU + uj 2 P(u). 

Let us for illustration consider some special cases. 
When the momentum is parallel to the magnetic field 
(6 = 0), the energies of the modes given in Eq. (29) 



reduce to |6|, u>e, and u> 2 /\b\, respectively. The type-II 
NG excitation is circularly polarized in the plane perpen- 
dicular to b, in the opposite direction than the gapped 
mode. This is in a nice analogy to the type-II NG boson 
in ferromagnets — the magnon — where the circular polar- 
ization corresponds to the Larmor precession of the elec- 
tron spin. Note that in the limit when uj t = ujg , the same 
conclusions hold regardless of the orientation of the mo- 
mentum. 

Second, when the momentum is perpendicular to the 
magnetic field (9 = 7r/2), the dispersion relations in 
Eq. (29) become \b\, uj t , and u>t<^e/\b\, respectively. The 



type-II NG mode is now polarized elliptically in the plane 
transverse to the magnetic field. The ratio of the princi- 
pal axes of the ellipse, oriented along and perpendicularly 
to the momentum, is given by uj t /ui(.. 

In the general already mentioned above, the 

polarization of the gapped mode is circular and trans- 
verse to the magnetic field, while that of the type-I NG 
mode is linear and aligned with the magnetic field, irre- 
spective to the orientation of the momentum. Only the 
polarization of the type-II NG excitation is affected by 
the momentum. 

In any case, there are only two (that is, d — 1) gapless 
phonons, and one of them has a quadratic dispersion re- 
lation. This can be explained by the combination of the 



conjecture ( 12 1 and the Nielsen-Chadha theorem. The 



non-commutativity of the components of the generator 
of magnetic translations suggests that there is one less 
phonon than one would naively expect. Then, to satisfy 
the Nielsen-Chadha theorem, at least one of the phonons 
must be type-II. 



C. Possible application 

The second-quantized version of the above model can 
be used to describe the rapidly rotating BEC of cold 
atoms. When viewed from the co-rotating frame, the 
Hamiltonian of the system picks an additional term, 
equal to — L ■ SI (coupling of the angular momentum L to 



the angular velocity SI), It then has formally the same 
form as the Hamiltonian of charged particles in the exter- 
nal magnetic field, provided the symmetric gauge is used 
and the replacement SI — > qB/(2m) made [31]. There- 
fore, the Hamiltonian of the rotating Bose gas has a two- 
dimensional continuous magnetic translational symme- 
try, in addition to the U(l) global symmetry correspond- 
ing to the conservation of particle number, as long as the 
trapping potential (modified by St) can be neglected. 

The formation of a triangular vortex lattice has been 
experimentally observed |32J, which implies spontaneous 
breaking of the continuous magnetic translational invari- 
ance. The collective mode due to the lattice distortion, 
called the Tkachcnko mode, has also been experimen- 
tally observed [33] . Similarly to our toy model, there 
is one gapless excitation with a quadratic dispersion re- 
lation and one gapped mode with the cyclotron gap, 
qB /m = 2SI [34] . The discussion in Section II E suggests 



there is still an independent Bogoliubov mode stemming 
from spontaneous breaking of the U(l) symmetry. We 
believe that our general arguments based solely on sym- 
metry properties may lead to a new understanding of NG 
modes in the vortex lattice. 



V. SUMMARY AND CONCLUSIONS 

In the present paper, we have developed a framework 
for the analysis of spontaneous symmetry breaking in the 
case that there is only a discrete translational invariance. 
While our prime motivation was to study spontaneous 
breaking of continuous translational invariance, our re- 
sults apply equally well to systems where this is broken 
explicitly by construction. 

We generalized existing theorems on NG bosons in 
quantum many-body systems, thus completing the pro- 
gram initiated in our recent paper [5] . Our main message 
is that as long as SSB is concerned, the usual division of 
symmetries into internal and spacetime ones is artificial. 
Instead, we propose that all symmetries that are uni- 
form in the sense defined in Section [IT] can, and should, 
be treated on the same footing. This naturally includes 
(uniform) internal symmetries as well as translational in- 
variance. 

To illustrate our general arguments, we analyzed a 
model for a (ferromagnetic) supcrsolid in two spatial 
dimensions using the mean-field approximation at zero 
temperature. Using numerical computations, we exam- 
ined the regime of weak coupling where the model fea- 
tures a homogeneous superfluid state, as well as the 
strong-coupling one where a supersolid crystalline order 
emerges. The properties of the spectrum of NG modes 
agree with our general conclusions in both phases. In par- 
ticular, we have demonstrated that the Bogoliubov mode 
and the longitudinal phonon in the supersolid state are 
distinct physical modes. 

Finally, we discussed the effect of an external magnetic 
field on the translational properties of the system. Using 
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a simple quantum-mechanical toy model, we argued that 
this changes the spectrum of NG modes qualitatively by 
making the dispersion relation of some of them quadratic 
at low momentum, and giving a gap to others. This issue 
would certainly deserve a more detailed investigation in 
a quantum-ficld-theoretic setting, with respect to its in- 
trinsic interest as well as the potential application to the 
physics of rotating Bose gases. 
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where we denoted q' = (g^fea). In the last step, we 
used the fact that, similarly to the coordinate decompo- 
sition X\ = X\ + R, any vector k\ from the reciprocal 
space can be uniquely expressed as a sum of a vector 
from FBZ and a vector G from the reciprocal lattice. 
As long as fci already lies in FBZ, the only solution to 
the condition qi — k\ — G = is qi = ki and G = 0. 
The second contribution to the commutator in Eq. ([3]), 
/ d d+1 x (0\(j) b (0)ja(x)\0) e i(ut ~ k x) , can be calculated in 
the same way and we therefore skip the details. 



Appendix A: Derivation of Eq. Q 

The derivation of Eq. ([3| is straightforward, if some- 
what tedious. It bases solely on the partition of unity 
translational property ([I]), definition of the averaged 
charge density in Eq. Q, and some elementary integrals. 
We would also like to first remind the reader of the Pois- 
son identity for the periodic (5-function in M. d , 

JW e *.K = E , d ' (fci _ G)j (A1) 
R G 

where G are vectors from the reciprocal lattice, G = 
y~li—i n{bi (rii € Z and a.; ■ bj = 2^6^). Taking now the 
first term in the commutator in Eq. pi), we obtain by a 
series of manipulations, 



e i(qi-ki)-R 



(A2) 
I 

Appendix B: Goldstone theorem in the effective 
action formalism 

In Appendix B of the paper [5] we provided a proof of 
the Goldstone theorem based on the quantum effective 
action in a form applicable to nonuniform symmetries. 
However, we for the sake of simplicity assumed contin- 
uous translational invariance of the ground state. Our 
aim here is to relax this constraint. Like in the rest of 
this paper, we will assume that the spontaneously broken 
symmetry is uniform; the extension of our argument to 
the fully general case is straightforward though. 

Consider a theory of a set of (not necessarily scalar) 
fields, <f>i{x), whose classical action is invariant under the 



d d+1 x (0\f a (x)M0)\0) 



i(ujt — k-x) 



d d q 

(2n) d 



(0\j a (t,x 1 ,x 2 )\£,q){Z,q\<f >b (0)\0)e 



i{ujt — h-x) 



[ dt [ d d ' Xl [ d d - d 'x 2 y\ [ 

J J J ^ J (FBZ 

/ dtJ2 I d d ' Xl J2 I ^% i (0\f a (t,x 1 +R,OM,q)^ q \M0)\0)e 
J ii uc £ Afbz) (27r) d 

/ dtJ2 I d d ' Xl J2 I ^r<0|i2(0,x 1 ,O)|e,*|')<e, 9 '|^(0)|0)e^- fa )-- B -*- ei 

J a Juc c JFBZ \ Zn ) 



i{ut — fel'SEl) 



d d - d ' X2 e 

—i[es(q')-u]t 



E / ^ 



FBZ 



d d ssi , 



no ^uc 



j° a (0, Xl ,0)e- ik ^ 



n 



(2 



'uc 
' R 



E / ^'gi(0|ja (fc) |e,q')(^g'l^(0)|0)27r ( 5( W -e 4 (g'))E <5d '(9i- fc i- G 
c Jfbz g 



= 2 7 r^(O| Ja(fc) |e,fc)(e,fc|0 6 (O)|O) ( 5( W -e ? (fc)), 
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d d+1 



V n tz^tttx FjlMv)] = 0, (Bl) 



infinitesimal transformation 5<fii(x) — 9Fi[(f>(x)]. Here 9 
is a parameter of the transformation and Fi is a local 
functional of the fields that does not depend explicitly 
on the coordinate x. (This is equivalent to the require- 
ment that the symmetry be uniform.) Provided that Fi 
is linear in the fields, the quantum effective action, T[<f>], 
shares the symmetry of the classical theory jTH] . We then 
obtain the invariance condition 

a 2 r[0 o ] 

where the integration measure involves time as well as 
d spatial coordinates and <poi(x) denotes the vacuum 
expectation value of the field in accord with the nota- 
tion introduced in Section UlI A 21 The second functional 
derivative of the effective action here represents the exact 
inverse propagator of the theory, (x, y). Assuming as 
in the rest of the paper a discrete translational invariance 
in R d and continuous one in W 1 ~ d as well as in time, this 
can be Fourier transformed as 

d d k 



+°° diu 

27T 



(FBZ) 



{2n) d 



(B2) 



x e 



i[ki 



c«--«,)+* a -(- a -»)]c-i( W> fc ; sB 1J y 1 : 



where fci lies in FBZ. Note that the Fourier transform 
of the propagator still depends on the coordinates within 
the Wigner-Seitz unit cell as a result of the lack of full 
continuous translational invariance. Thanks to the as- 
sum ed p eriodicity of the ground state, <poi(y) — < Poi(yi)t 
Eq. (Bl ) becomes 



d d m 

r. 



grMo,0-x 1 ,y 1 )F j [My 1 )] = 0. (B3) 



For a translationally invariant ground state, Fj[<fio(y~i)] 
is a constant and it follows immediately that the inverse 
propagator must have a gapless pole at zero momen- 
tum once the symmetry is spontaneously broken, that 
is, -F}[</>o] is nonzero. However, in presence of a spatially 
varying order parameter, the implications of Eq. (B3) for 
the excitation spectrum are less straightforward. 

It is useful to first recall the spectral representation of 
the propagator. Using the definition of the propagator in 
the operator formalism, iQij(x,y) = (0\T{(f>i(x)(j)j(y)}\0) , 
and inserting the partition of unity ([2| , one arrives at the 
following expression for the Fourier transformed propa- 
gator, defined analogously to Eq. ( |B2| , 

E 



Gij(u,k;xi,yi) 

"(0|^(Q,xi,0)|e,fe)^,fe|^(0,yi,0)|0) 



(0\4>j(0, yi, 0)|£, -k)(£, -k\<f>i(P, x u 0)\0) 



(B4) 



-k) - i0 



This ensures that even with just a discrete transla- 
tional invariance, one can still extract the quasipar- 
ticle spectrum from the poles of a conveniently de- 
fined propagator. To conclude our argument, we note 



that the condition of Q 1 being inverse to Q reads 
/ d d+1 zg ik (x,z)G^ J 1 (z,y) = Sij5 d+1 (x - y), which takes 
the following form in the frequency-momentum space, 



d d Z! 



g ik (uj, k; xx,zi) g k Hu, k; Zi, yi) 



n ^ „ KV -,.-,- ± ,- ± ,* kj (B5) 

= <5ijfi uc (5 d (xi - yi). 

Multiplying this equation by Fj[(f> (yi)] a nd in tegrating 
over yi, we assert with the help of Eq. (B3) and the 



spectral representation (B4) that as soon as Fj[(f>o(xi)] 



is nonzero at some x±, the propagator (at the same point) 
must have a pole such that £g(0) = 0. 

Finally, let us remark that the proof would have been 
technically much simpler had we considered the effective 
action (and in turn the propagator) as a functional of 
fields spatially averaged over the Wigner-Seitz unit cell, 



i(t,R, x 2 ) = 



d d xi 



4> i (t,x 1 + R,x 2 ). 



(B6) 



The order parameter <j> 0i would then be completely time 
and coordinate independent and the left-hand sides of 
Eqs. (B3| and (B5| would become mere products in mo- 



mentum space. The price for this simplification would be 
that we could make no conclusions about the spectrum in 
the cases that the symmetry is spontaneously broken but 
the spatial average of the order parameter vanishes. The 
above given proof is more general since it only assumes 
that the order parameter is nonzero at some point. 



Appendix C: Calculational methods 

In this appendix, we collect some details of the tech- 
niques used in the numerical computations in Section [HT 

1. Minimization of the energy 



To minimize the energy functional (16), we exploit 



the assumed periodicity of the classical field and per- 
form the Fourier decomposition of the field, ipc = 
J nc d d x ?p(x)e~ lG ' x . The total energy divided by the 
space volume is then equal to the energy density averaged 
over the unit cell, which is given by 



5 E ^iG^iG-qViqWjc^jG' 



(CI) 



+<?■ 



G.G'q 



Thanks to the parity invariance of the Lagrangian, which 
is naturally assumed to be shared by the ground state, we 
can restrict our attention to even-parity classical fields, 
ipo(x) = ipo(-x). This, together with the reality of the 
classical field, implies the relation ipc = ip—G — V'g' 
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which significantly reduces the number of variables in 
the minimization procedure. 

The minimization of the averaged energy density is 
accomplished by first fixing the primitive lattice vec- 
tors {a,i}f =1 and minimizing £[ip] with respect to the 
Fourier components {V'g}- This determines the function 
£({a,i\) = £[-0o({cti})], which is subsequently minimized 
with respect to {a{\. The actual crystal form of the clas- 
sical field ^o(x) is thus found. 



where the (infinite-dimensional) matrices K, T and W 
are defined by 



K. 



(fc + Gi) 
2<7 



^GiG 2 + U G 1 -G 2 i 



Tg 1 g 2 = iu(k)5 GlG2 , (C3) 

W Gl G 2 = ^2 4>0G 1 -G>i>0G>-G 2 v(k + G'), 



G' 



Calculation of the band structures 



The Schrodinger equation (21) can be rewritten in 
Fourier space as 



T G 1 G 2 ( PiG 2 = -{K Gl G 2 +2d i xWG 1 G 2 )<t>iG2, 
TG 1 G 2 (f>iG 2 = —Kg 1 G 2 ( PiG 2 , 



(C2) 



and ug = v{G) ^2 G > ipoG-G'^OG 1 - Because Eq. (C2) 



represents a set of homogeneous linear equations for the 
Fourier components 4hg and <piG (for each fixed i), the 
band structure can be found from the condition that the 
determinant of the matrix of coefficients in these equa- 
tions be zero. 
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